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Abstract 



! We discuss the time-dependent CP asymmetry of decay B — ► in an extension of the Stan- 



dard Model with both two Higgs doublets and additional fourth-generation quarks. We show that 
although the Standard Model with two-Higgs-doublet and the Standard model with fourth genera- 
tion quarks alone are not likely to largely change the effective sin 2(5 from the decay of B — * 4>Ks, 
the model with both additional Higgs doublet and fourth-generation quarks can easily account for 
the possible large negative value of sin 2/? without conflicting with other experimental constraints. 
| In this model, additional large CP violating effects may arise from the flavor changing Yukawa in- 

teractions between neutral Higgs bosons and the heavy fourth generation down type quark, which 
can modify the QCD penguin contributions. With the constraints obtained from b — * sss processes 



such as B — > X s j and Amgo, this model can lead to the effective sin 2(3 to be as large as —0.4 in 
\ the CP asymmetry of B — > 4>Ks- 
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I. INTRODUCTION 



With the successful running of two B factories in KEK and SLAC, precise measurements 
of the time-dependent CP asymmetries as well as the directly CP asymmetries in rare B 
decays become available. Among those interesting decay modes, the most important one, the 
CP asymmetry of B — > J / ipKs has been successfully measured, and a very good agreement 
with the Standard Model (SM) prediction on sin 2(3 was found. 

However, the recent Belle results on sin 2(3 from B — > <(>Ks, although with significant er- 
rors, have indicated that the value of sin 2(3 from different decay modes could be significantly 
different. The most recent measurements give 



sin 2(3 
sin 2(3 



0.47 ± 0.341^ (Babar), 
-0.96±0.5l°;°i( B elle). 



(1) 



Of course, it is too early to draw any robust conclusion from the current preliminary data. 
Nevertheless, it opens a possibility that large new physics effects may show up in the b — > sss 
processes, which has already triggered a large amount of theoretical efforts in examining the 
possible new physics contributions from various models. Besides the models related to 
supersymmetry which are the most promising ones, there are also a large class of models 
based on simple extensions of the matter contents of the SM, such as the standard models 



with two-Higgs-doublet (S2HDM) 0, 13, II B 0, B H E EUii.HJliJll 



19, 17, UM and 



the standard model with fourth-generation fermions (SM4)|l9l. I20LI21. 2jJ|23j etc. However, 
the most recent studies have pointed out that the contributions from the above mentioned 
two types of models to B — > <pKs are in general not large enough to account for a large 
negative value of sin 2(3 in B — > <pKg ( for example sin 2(3 ~ —0.5 ) Eil liH |26| . 

In this paper, we show that although due to the constraints from other experiments 
such as 6 — * S7 and Am# etc., the general S2HDM and the SM4 alone are not likely to 
largely change the effective sin 2(3 in B — > 4>Kg, a model with both an additional Higgs 
doublet and 4th-generation quarks (denoted by S2HDM4) can significantly change the value 
of sin 2(3 without contradicting with other experimental constraints. In this model, new large 
CP violating contributions may arise from the flavor- changing Yukawa interactions between 
the neutral Higgs boson and the 4th-generation down type quark b' (with my ^> nit,), which 
changes the Wilson coefficients for QCD penguin operators and results in a large modification 
of effective sin 2(3. This mechanism is different from the case in the S2HDM in which the 
dominant contribution comes from changing the Wilson coefficients of the elect ro(chromo) - 
magnetic operators. The latter is subjected to a rather strong constraint from b — > 57 and 
therefore can not give enough contributions. 

Let us begin with some model independent discussions. The definition of effective sin 2(3 
in B — > (pKs is 



sin 2(3, 



elf 



Im 




= Im 




A 





1 +re- ie ) 



A su (l + re+ ie ) 



(2) 



where (3 is the SM value with sin 2(3 = 0.7151^ |27|. -4 S m(*4 S m) is the SM value of the 
decay amplitude of 5° (B°) — > 4>Ks- Here two parameters r and 6 parameterize the relative 
size and the additional CP violating phase of the new physics contributions. To get an idea 
of how sin 2/5 e fj is changes with the new physics contribution, we take some typical values of 
the phase 9, calculate the values of sm2(3 e s, and shown them in FigJT] 
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FIG. 1: The value of sin2/3 c fj as a function of r. The solid, dashed and dotted curves corresponds 
to 9 = 7r/2,7r/3 and 7r/6 respectively. 



As it is shown in the figure, to explain the possibly large negative sin2/3 e ff, for instance, 
close to —0.5, in the case that 9 is maximum (7r/2), the value of r should be close to unity. 
For smaller 9 such as 7r/3 and 7r/6, the value of r must be even larger. Therefore, to generate 
a large negative value of sm.2(3 e ff in the range of —0.5 ~ —1.0, the magnitude of the new 
physics contributions must be as the same order of magnitude as the one in the SM. 

However, the new physics contributions must be constrained by other experiments, es- 
pecially by the b —>■ s transition related processes. The most strict constraint comes from 
the radiative decay 
(3.28±j}:£) x io- 4 E 



of B ->• X.<y. The current data of Br(S -> X sT , £ 7 > 1.6 GeV) = 
I 2^, 3(J EH is well reproduced in the frame work of the next-to-leading 
order calculations in the SM (see e.g. |32|, |33[). Thus, if the new physics contribution carries 
no new phase, there is very little room for the new physics parameters. But in the case 
that new phases present, the parameter space could be enlarged. This is because the data 
of B — > X s 7 only constraints the absolute values of the Wilson coefficient CV 7 , if the new 
physics contribution does not change the absolute value of C 77 , there will not be a serious 
problem. Thus the following relation must be satisfied for any new physics model 



G 



77 1 



I^-ySM i ^NEWI 
IW7 "+* W7 



°7 7 I 



(3) 
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with Cf^ 1 and C^ EW being the effective Wilson coefficient evaluated at the low energy scale 
(/i « rrib) from SM and new physics models respectively. In this case, the absolute value of 
C^L could vary largely from close to zero to about —2C^r, which seems large enough for 
explaining the CP asymmetry in B — ► 4>Ks- However, it follows from Eq.fJHJ) that the data 
onB-> X s 7 do strongly constrain the form of C^ y EW , namely, the new physics must interfere 
in such a way that the total effect is roughly equivalent to adding a phase factor to Cf^ 1 , i.e 
C*7 7 ~ | Cf^ 1 1 e* . Let us take an illustrative example in which the new physics contribution 
is purely electro (chromo)-magnetic and satisfy C 7l = ICj^le 10 and also C 8g = \C^\e td at 
the scale of mjy. Varying the value of 9 from to 2tc and then running down to the low 
energy scale of \x ~ through renormalization group equation, one finds that the value of 
sin2/9 e fj in decay B — > <pKs only changes from 0.5 to 0.8, This naive discussion shows that if 
the dominant contribution from a new physics model is coming from C^ff^, the change to 
sin2/5 e fr from the its SM value is limited. Unfortunately, the S2HDM belongs to this class 
of model. The recent analysis have confirmed that within S2HDM, the value of sin2/3 e ff can 
reach zero, but not likely to be largely negative jiil. liH l26l| . 

For the model of SM4, there are additional up {if) and down (&') type quarks. The new 
phases may come from the extended Cabbibo-Kobayashi-Maskawa(CKM) matrix which is a 
four by four matrix in this model and contains undetermined matrix elements of V t , , V q y etc. 
To avoid the precise data of electro-weak processes, the mass of b' {f ) has to be pushed to 
greater than ~200 GeV( ~ 300 GeV). However, phenomenological study showed that with 
the constraint of B — > X s 7 and B° s — B° s mixings being considered, its contribution to the 
CP violation of B — > 4>Ks is not large enough either [2]]]. Thus if the large negative value of 
sa\2f3 e ff in decay B —> <pKs is confirmed by the future experiments, the above mentioned 
two models ( i.e. S2HDM and SM4 ) will not be favored. 

sectionThe model of S2HDM4 

There are several directions in constructing models beyond the SM, such as enlarging the 
gauge groups to SU(5), SU(10) andEg etc., introducing new symmetries like various SUSY 
models, and expanding the matter contents, i.e., more fermions and Higgs bosons. The 
models of the last type can be regarded as simple extensions of the SM which keep the same 
gauge structure but still have rich sources of new contributions. The typical ones are the 
above mentioned S2HDM and SM4. 

In this paper we would like to a step further to consider a model with both two-Higgs- 
doublet and fourth-generation quarks (S2HDM4). In this model, there are new Yukawa 
interactions between Higgs bosons and heavy fourth-generations quarks. Since in general 
the Yukawa interaction is expected to be proportional to the coupled quark mass, the new 
Yukawa couplings are much stronger than that in the S2HDM and SM4 . Unlike in the case 
of S2HDM, where the b quark contribution to the QCD penguin diagram through neutral 
Higgs boson loop is strongly suppressed by the small b quark mass, the same diagram with 
intermediate b' quark may significantly contribute to the related processes |34|. This new 
feature only exists in this combined model, and is of particular interest in studying the CP 
violation of B — > <pKs and other penguin dominant processes. 

The Lagrangian for the S2HDM4 is given by 

Cy = jiXFfauR + $ L Yf<f>xd R + $ L YVfau R + ^ L Y 2 D <j)2d R + H.c (4) 

with the extended quark content of u^r = {u, c, t, t')L,R an d ^l,_r = {d, s, b, b')^R. The 
Yukawa coupling matrices are 4-dimensional matrices accordingly. The two Higgs 

fields 0i,02 have vacuum expectation values (VEV) of vie tSl and f 2 e l<52 respectively, with 
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a/I^iI 2 + |fi| 2 — v — 246 GeV . The relative phase S = 61—62 between two VEVs is physical 
and provides a new source of CP violation 0, 0, ■ In the mass eigenstates, the three 
physical Higgs bosons are denoted by H°, A , and H ± respectively. Due to the non-zero 
phase 6, all the Yukawa couplings become complex numbers in the physical mass basis, even 
they are all real in the flavor basis. For simplicity, throughout this paper, we assume that 
the CKM matrix elements associating with t', i.e. V t > q are ignorablly small and will only 
focus on the neural Higgs boson contributions. 

In the mass basis, the Yukawa interactions between neutral Higgs bosons and quarks have 
the following general form 

£y = VijQiLqjR(p + H.c, (5) 
with = H° or A . The Yukawa coupling rjZ is usually parameterized as 



./m ai m a . 

4 = (6) 

In the Chen-Sher ansartz [35| motivated by a Fritzsch type of Yukawa coupling matrix, the 
values of all £ 3i? s are of the same order of magnitude. However, from other textures of the 
coupling matrix the relations among £g <( ^s are different 113, HH]. In the general case, they 



should be taken as free parameters to be determined or constrained by the experiments. 
The effective Hamiltonian for AB = 1 charmless B decays reads 



H cS = % [K h V^(C?Ql + C^Ql) + V ch V*(C1Ql + C C 2 Q 



V2 



10 



,i=3 



(7) 



where the operator basis QiS can be found in Ref . |39f| . In this model, the relevant Wilson 
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coefficients at the scale of mw from this model is given by 

Ua,(M w ) 



Ci(M, 



w i 



2 4tt ' 

11 a s (M w ) 35 a e 



W)=l- 6 4I 18 4.' 

c 3 (m w) = -f^t^t.o+i^d^TO + ^^a.^^M), 

+ ^-^— (2£ (* t ) + Co(^)), 

07T Sin #14/ 

rfM \ ®s(M w ) - 2 m myy/mbms m 

C 4 {M W ) = — — — (£ (a*) + |6t| £ (l/) + 0T/ T ~ 2 £wtsb>E (y)), 

C 5 (M W ) 247T ' ^ 21/ y/m 2 & 

C 6 (M W ) - — C-ooW + lcttl ^0 UO + ot/ t/*^,2 wUAi 

8^ ^bW™t 

^ x A{x t ) 1 ( 2 my^/m b m s 



+B(y)\^ b \e -B{y) g— ^U', 

z ' v tb v ts" l t" L b 

+(z/)l&6|e - E (y) 9 ir v . - ZvbLv, (8) 

with a s (mw) and a em being the strong and electro-magnetic couplings at scale mjy. The 
mass ratios a; t , 3/ and y' are defined as x t = m^/m^y, y = m^/m 2 H± and y' = Wy/ mL 
respectively. The loop integration functions are standard and can be found in Refs. p, |4(j, 
4l|,|43. Here we have ignored the coefficients for the electro-weak penguin diagrams since 



their effects are less significant in the decay of B — > <pKs- 

Note that the new contributions to QCD and electro (chromo)-magnetic operators de- 
pends on different parameter sets. In the QCD penguin sector, the contribution depends on 
Qb'isv where in electro (chromo)-magnetic sector it depends on both i b i b isb' and Qyisv- It is 
convenient to define two weak phases 9\ and 62 with 

CwU* = \Zbb>U¥ ei and UU> = \UU>\e ld2 - (9) 

Since in general and Q b , are complex numbers and £, b , b 7^ Q b ,, the two phases are not nec- 
essary to be equivalent. The presence of two rather than one independent phases is particular 
for this model, which gives different contributions to the QCD penguin and electro(chromo)- 
magnetic Wilson coefficients. The interference between them enlarges the allowed param- 
eter space. Note that the Wilson coefficient for QCD penguins may be complex numbers 
which provides additional sources of CP violation. To make a comparison, let us denote 
the Wilson coefficients in the SM by Cf M . Taking £ Vb = £ sV = 0.8, 9 X = 0.5, 6 2 = -1.2 
and ra#o = my = 200GeV as an example, in the range of 40 < < 60, the ratio of 
Cz/C^ M (C A /Cf M ) has an imaginary part between -0.27 and -0.4(-0.6 and -0.8). These large 
imaginary parts plays an important role in CP violation. 



6 



II. CONSTRAINTS FROM B -» X s7 ANDB^ - 5° MIXING 

Before making any predictions, one first needs to know how the new parameters in this 
model are constrained by other experiments. For the process we are concerning, the most 
strict constraints comes from b — > sss processes such as B — > X s 7 and B® — B® mixing, etc. 

The expression for B — > X s 7 normalized to B — > X c eu e reads 



Br(E -> X sl ) 
Br (5 -> X r eu P ) 



6|VtbVS 


2 a 


7T 


Kb 


2 f{m c /m b ) 



CrM\ 2 (10) 



with/(z) = l-82 2 -242 4 ln^+8z 6 -z 8 and Br(5 -> X c ez/ e ) = 10.45%. The low energy scale 
\i is set to be m&. Using the Wilson coefficients at the scale mjy and running down to the m& 
scale through re- normalization group equations, we obtain the predictions for Br [B — > X s j). 
For simplicity, we focus on the case in which the h' contribution dominates through H° loop, 
namely, we push the masses of the charged Higgs if ± and the other pseudo-scalar boson 
A to be very high (m#± , m^o > 500 GeV) and ignore their contributions. We take the 
following typical values of the couplings 

1 6b' I = 50, = 0.8, = 0-8, and m H o =m' h = 200 GeV, (11) 

and give in FigOthe value of Br(B — > Xg^y) as a function of 9\ with different values of 62- 

From the figure, one finds that two separated ranges for parameters Q\ and 82 are allowed 
by the data 

-1.4 < 6 2 < -1.2 and 0.4 < 6 2 < 0.7 for 0.5 < Q x < 1.5, (12) 

Note that we do not make a scan for the full parameter space, nevertheless the above 
obtained range are already enough for our purpose. Among the two allowed ranges, the one 
with —1.4 < 82 < —1.2 is of particular interest. It will be seen below that in this range, 
the contribution to the CP asymmetry in B — > <pKs could be significant. In FigOU we also 
give the allowed range of Q\ with difference values of 62- One finds that the allowed range 
for B\ is larger compared with 6 2 . In this figure, the interference between two phases 9\ 
and 62 is manifest. For 8 2 in the range of (—1.0, —0.8), the allowed value for 9\ is a narrow 
window around zero. But for 9 2 in the range of (—1.4, —1.2), the allowed range for B\ could 
be between 0.5 and 2.0. Compared with the S2HDM in which only one phase appears, 
this interference effect for two phases enlarges the parameter space under the constraint of 
B — > Xg'y. Thus large contributions to the other processes is possible in this model. 

The other b — > sss process which could impose strong constraint is the mass difference of 
neutral B° s meson. The measurements from LEP give a lower bound of Am# s > 14.9ps _1 . 
In this model, the b' contributes to AmRonlv through box-diagrams. The box diagram 
contribution to Am# s is given by BQMIII 



bn z 4 



+ 2 Vu i/tlsttl B v u/t> Vw) + -rVtt y { 0T/ T7i 2 wwJ B v yy)s \ 16 > 

1 A V tb V ts m t 



where Gp = 1.16 x 10 5 GeV 2 is the Fermi constant. fs a and B Bs are the decay constant 
and bag parameter for In the numerical calculations, we take the value of fs s \/B Ba = 



7 




FIG. 2: The branching ratio of B — > X s ~f as a functions of #2 in the model of S2HDM4 . The solid, 
dashed and dotted curves correspond to 9\ = 1.5, 1.0 and 0.5 respectively. Other parameters are 
taken from Ea. (|lll) . 



0.23GeV. r/ijS are the QCD correction factors. The loop integration functions of B, u u !m 



(V) 



can be found in R,efs|4lL l42t l43{ |. The mass ratios are denned as y t = mf/m 2 H± ,y w = 
m t/ m w an d y' — m b>/ m 'H0 respectively. Note that in the mass difference of B® mesons, the 
contribution from S2HDM4 only depends on the parameter So, only the phase Q\ 

will present in the expression. 

Using the above obtained typical parameters in Eq. ()ll|) . the contribution to Ame s is 
calculated and plotted as a function of Q\ in FigEJ The figure shows that the current data 
of ArriB s do not impose strong constraint on the value of Q\. 

The neutron electric dipole moment (EDM) is expected to give strong constraints on the 
new physics. In the SM, the neutron EDM is zero at even two loop level. The current experi- 
mental upper limit gives EDM< 1.1 x lO _25 ecm.0]. In general, the new physics contributes 
to the neutron EDM through one loop diagrams. In the presence of new scalars, additional 
significant contributions may arise, for example from the Weinberg gluonic operator 45] and 
also the two-loop Barr-Zee type diagrams jia, etc. 



However, we note that all the above three type of mechenisms are not related to b — > s 




FIG. 3: The branching ratio of B — > X s ~f as a functions of 9\ in the model of S2HDM4 . The solid, 
dashed, dotted and dot-dashed curves correspond to 62 = 1.4, 1.2, 1.0 and 0.8 respectively. Other 
parameters are taken from Ea. ()lll) 



flavor- changing transitions and therefore will involve different parameters in this model. 
For the one-loop diagrams, the neutral EDM is mostly related to and £ t {b') through 
u{d)— quark EDM. For Weinberg three gluonic operator, the dominate contribution is from 
intereral b' loop. Thus it is related to Similarly, for two-loop Barr-Zee diagram, the 
b'— quark loop will play the most important role and the couplings involve only £, u (d),t,b'b> 
etc. Thus the neutron EDM will impose strong constraints on other paramerters in this 
model and has less significance in current studying of decay B — > <pKs- This is significantly 
different from the S2HDM case in which the t— quark alway domains the loop contribution 
and the couplings £ tt an d £&& are subjected to a strong constraint from neutron EDM. 

Other constraints may come from K° — K° and — B® mixings. But those processes 
contain additional free parameters such as the the Yukawa coupling of and £ s &/, the 
constraints from those processes are much weaker. 



9 




-3 -2-10 1 2 3 

FIG. 4: The B® meson mass difference AB® as a function of 9\ in the model of S2HDM4 . Other 
parameters are taken from Eq.(|ll|). The shadowed region is excluded by the data of AB®. 

III. CP ASYMMETRY IN B -> cpK s 

Now we are in the position to discuss CP asymmetry in B — > <pKs- The decay amplitude 
for B — > (f)K°rea.ds 

G 1 
A(B° d -> c/)K ) = --I=V*V th (a 3 + a A + a 5 - ~{a 7 + a 9 + a 10 ))X, (14) 

with X being a factor related to the hadronic matrix elements. In the naive factorization 
approach X = 2/^m^(e • pB)Fi(m<f>), where e, pb, F\ are the polarization vector of 0, the 
momentum of B meson and form factor respectively. The coefficients a« are defined through 
the effective Wilson coefficients Cf ff s as follows 

Since the heavy particles such as H ±,Q , A and b' has been integrated out below the scale of 
mw-, the procedures to obtain the effective Wilson coefficients are exactly the same as 
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in SM and can be found in Ref.[4£ 

Using the above obtained parameters allowed by the current data, the prediction for the 
time dependent CP asymmetry for B — > <pKs are shown in Fig|S] 




FIG. 5: The prediction for sin 2f3 e s as a functions of 0\ with different value of 62- The solid, dashed, 
dotted and dot-dashed curves corresponds to 62 = —1-4, —1.2, —1.0, —0.8 respectively. 

In the figure, we give the value of sin2/3 e ff as a function of 9\ with different values of 
#2=1-4,1.2,1.0 and 0.8. Comparing with the constraints obtained from B — > X s 7 and B®—B® 
mixings, one sees that in the allowed range of —1.4 < 82 < —1.2 and 0.5 < Q\ < 1.5, the 
predicted sin 2/3 e fr can reach —0.4. 

It is evident that the large negative value of sin2/? e jj is a consequence of the interference 
effects between Q\ and 9 2 and therefore is particular for this model. For zero value of 61, 
there is no new phase in the QCD penguin sector. From FigJHJ the allowed range for 9 2 is 
— 1.0 < 6 2 < —0.8. Then, it follows from FigEl that in this range the predicted sin2/3 e jj is 
at around zero. But for 9% ~ 0.5, the allowed range for 82 is changed into —1.4 < 8 2 < —1.2 
and the predictions for sin2/3 e // is much lower in the range of (—0.4, —0.25). 
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IV. CONCLUSIONS 



In conclusion, we have discussed the CP asymmetry of decay B — > (j)Ks, in the model of 
S2HDM4 which contains both an additional Higgs doublet and fourth generation quarks. In 
this model, since the fourth generation b' quark is much heavier that b quark, the Yukawa 
interactions between neutral Higgs boson and b' is greatly enhanced. This results in signifi- 
cant modification to the QCD penguin diagrams. We have obtained the allowed range of the 
parameters from the process of B — > X s ^ and ArriB a - Due to the more complicated phase 
effects, in this model the constraints from those process are weaker than that in S2HDM 
and SM4. The effective sin2/3 e ff in the decay B — > (pKs is predicted with the constrained 
parameters. We have found that this model can easily account for the possible large negative 
value of sin 2f3 without conflicting with other experimental constraints. 

In this paper we focus on the case in which H° domains. It is straight forward to find 
that the contribution from the other pseudo-scalar A follows the same pattern. In the 
case of small mixing among the neutral scalars, the Yukawa couplings for H° and A are 
directly related [l2l ]. We find that for m A o m 200 GeV <C m H o its contribution to the decay 
amplitude of B — > <pKs is similar to the case of the H° dominance discussed above. For 
the case that m^o is close to m#o, the contribution from them are comparable, and the 
interference between the two could be important. 

Since this model contributes new phases to QCD penguin diagrams, it remains to be 
seen if it has sizable effects on other penguin dominant processes, such as in the hadronic 
charmless B decays. Similarly, it is expected that in this model there are also significant 
contributions to the electro-weak (EW) penguin diagrams which deserves a further inves- 
tigation (for recent discussions on EW penguin effects on B — ► <f)K see, e.gji^, IHlf.) 
It is well known that the EW penguin plays important roles in rare B decays. The cur- 
rent data on B — > -7T7T, ttK have indicated some deviations from results based on the SM 
53L 54l 56l 57^ . It is of interest to further investigate the new physics contributions 



to those decay modes within this model. 
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